#!/usr/bin/env python

from math import pi, e, sin, cos, sinh, cosh

def exact_solution(x,y):
	return y*(1-y)*x**3


